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A variational model is proposed to describe the magnetic 
properties of type-II superconductors in the entire field range 
between Hci and IIc2 for any values of the Ginzburg-Landau 
parameter k > l/\/2. The hexagonal unit cell of the trian- 
gular flux-line lattice is replaced by a circle of the same area, 
and the periodic solutions to the Ginzburg-Landau equations 
within this cell are approximated by rotationally symmetric 
solutions. The Ginzburg-Landau equations are solved by a 
trial function for the order parameter. The calculated spatial 
distributions of the order parameter and the magnetic field 
are compared with the corresponding distributions obtained 
by numerical solution of the Ginzburg-Landau equations. The 
comparison reveals good agreement with an accuracy of a few 
percent for all k values exceeding k « 1. The model can be 
extended to anisotropic superconductors when the vortices 
are directed along one of the principal axes. The reversible 
magnetization curve is calculated and an analytical formula 
for the magnetization is proposed. At low fields, the theory 
reduces to the London approach at k ^ 1, provided that the 
exact value of Hd is used. At high fields, our model repro- 
duces the main features of the well-known Abrikosov theory. 
The magnetic field dependences of the reversible magnetiza- 
tion found numerically and by our variational method prac- 
tically coincide. The model also refines the limits of some 
approximations which have been widely used. The calculated 
magnetization curves are in a good agreement with experi- 
mental data on high-Tc superconductors. 

PACS numbers: 74.20.De, 74.60.Ec, 74.25. Ha 



I. INTRODUCTION 

The solution to the Ginzburg-Landau (GL) equations 
found by Abrikosov was used widely |2||] to study the 
properties of type-II superconductors at low and high ap- 
pUed magnetic fields H , i.e. close to the lower and upper 
critical fields Hd and Hc2 , respectively. At H ^ Hd the 
intervortex spacing is much larger than the vortex core 
size if the GL parameter is large, k ^ 1. Therefore in the 
London model, which is commonly used at low fields [Q , 
the order parameter in the superconductor is assumed to 



be constant. In this case, the flux- line lattice (FLL) can 
be treated as a set of independent vortices: the magnetic 
flux density is a linear superposition of the contributions 
of individual vortices, and each contribution coincides 
with the field of an isolated vortex. The energy of the 
system is the sum of the self-energies of the vortices and 
their pairwise interaction ||^. At high fields, the Lon- 
don model looses its applicability because the fraction of 
the total volume of superconductor occupied by the vor- 
tex cores is no longer small ||^ . Different approximations 
were proposed to include the vortex cores and to extend 
the applicability limits of the London model |^ . 

The problem of solving the GL equations for the ide- 
ally periodic FLL can be simplifled considerably by re- 
placing the hexagonal unit cell of the vortex lattice by 
a circle of equal area (Wigner-Seitz approximation). In 
this approach both the order parameter and the mag- 
netic flux density within the cell have axial symmetry. 
The presence of other vortices is taken into account by 
the boundary condition: the supercurrent density equals 
zero at the cell boundary. This method was used in Ref. 
[Q , where an explicit expression for the magnetization in 
low fields was found, and the results are in good agree- 
ment with the predictions of the London model ||^,^ . In 
Refs. [p[-jl3|, the GL equations and the equations of the 
microscopic theory of superconductivity were solved nu- 
merically in the framework of circular cell approaches 
and it was shown that this approximation not only yields 
good results at low induction but also a.t H ^ Hc2- A 
similar approach was used in Ref. [2l| ] , where the system 
of GL equations was reduced to a single equation which 
can be solved numerically. The Wigner-Seitz approxima- 
tion can be extended even to the case of exotic pairing 
symmetries, when the order parameter has two complex 
components [Q. 

A numerical method to find the periodic solutions to 
the GL equations was developed in Refs. |l^,|l^. This 
exact method accounts for the actual symmetry of the 
vortex lattice. It allows calculating the spatial distri- 
butions of the magnetic field and the order parameter 
within the unit cell, the elastic shear modulus of the FLL, 
and the magnetization for any FLL symmetry and any 
induction B and GL parameter k with any desired ac- 
curacy. However, until recently there was no adequate 
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approach allowing one to find the magnetization analyt- 
ically in the entire field range Hd < H < Hc2 and to 
obtain explicit formulas which may be used to analyze 
experimental data. 

In Ref. Clem proposed a model to solve the GL 
equations using a trial function for the order parameter 
(or GL function) \ip\: 



\^\ = f{r) 



(1) 



where S,v is a variational parameter. This model yields 
an approximate explicit expression for the local magnetic 
field of an isolated vortex: 



hir) 
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(2) 



where Kn are modified Bessel functions. Here and be- 
low the following dimensionless variables are used: dis- 
tance r, magnetic flux density h, and order parameter / 
are measured in units of A, Hc\/2, y^—a/ (3, respectively, 
where A is the London penetration depth, He is the ther- 
modynamic critical field, and a and /3 are the GL coef- 
ficients. In this notation, we have Hc2 = — 27t/k, 
where $o is the magnetic flux quantum. For k 3> 1 the 
minimization of the free energy gives ^i, ~ \/2/k and 



Hci = — (luK + s), 



(3) 



where e — 0.52. The exact value e ~ 0.50 was calculated 
in Ref. ||l8| (see also Ref. from the numerical solu- 
tion to the GL equations for an isolated vortex. Note that 
the lower critical field cannot be found self-consistently 
in the framework of the London model, because in this 
approach the magnetic flux density diverges on the vor- 
tex axis. Therefore, Hd should be regarded as a free 
parameter in the London expression for the magnetiza- 
tion However, the approximation for Hd found by 
cutting off the field of a vortex at a distance equal to the 
coherence length ^, is often used in the London approach 
as well ||2^. Note that recently Clem's trial function (1) 
was applied to the study of the vortex core structure in 
superconductors with mixed {d~\-s) two-component order 
parameter |2^]. 

Hao et al. |23|| (see also Ref. ^J]) extended the model 
pTf to larger magnetic fields up to Hc2 through the linear 
superposition of the field profiles of individual vortices. 
In this model, the trial function (1) is multiplied by a sec- 
ond variational parameter f^o to account for the suppres- 
sion of the order parameter due to the overlapping vortex 
cores. This model enabled the authors to calculate 
the magnetization of type-II superconductors in the full 
range Hd < H < Hc2- Their analytical formula is in a 
good agreement with the well-known Abrikosov high-field 
result. For the case of low fields, Hao and Clem argued 
psf that the London model is quantitatively incorrect 
(it does not give the correct asymptotics at iJ ^ Hd) 
since the contribution of the vortex cores to the total 



free energy could not be taken into account in this ap- 
proach. The Clem-Hao model was further extended to 
include anisotropy p^ , |2^ . This approximation is now 
widely used for the analysis of the experimental data on 
magnetization of type II superconductors p^-| 2^. 

However, it has been recently shown in RcfT pCfl that 
the Clem-Hao model has some drawbacks. It was argued 
that the procedure of obtaining the local magnetic flux 
density by a linear superposition of contributions of indi- 
vidual vortices in the form used in Ref. is valid only at 
low fields. The application of this approach to the entire 
field range Hd < H < Hc2 leads to an appreciable dis- 
agreement between the Clem-Hao model and Abrikosov's 
high-field result |^ . In the original papers of Hao et al. 
p3| , p5|j26| , this disagreement was made up by the use of 
a non-selfconsistent field dependence of the variational 
parameters. Then, in calculating the magnetic free en- 
ergy in Ref. [ p3[ the lattice sum was approximated by 
an integral. As it was shown in Ref. |^ , this procedure 
leads to a noticeable error in the magnetization at low 
fields. This error and the use of an inaccurate value of 
Hd have led the authors |^,^ to the conclusion about 
the quantitative incorrectness of the London approach at 
low fields. Note that in Ref. |2|] it was argued too that 
the Clem-Hao model overestimates the effect of suppres- 
sion of the order parameter in the vortex cores at small 
fields. 

In this paper, we propose a variational model for the 
description of the regular fiux-line lattice in a more con- 
sistent fashion as compared to Ref. . Our variational 
procedure is based on Clem's trial function (1). How- 
ever, in contrast to the Clem-Hao model, we do not use 
the superposition of vortex fields. Instead, we apply the 
circular cell method and calculate the magnetic flux den- 
sity directly from the second GL equation. The model 
enables us to find analytical expressions for the local mag- 
netic field and the order parameter. The results of our 
variational procedure are compared with the results of 
the numerical solution of the GL equations. This com- 
parison reveals that the analytical formulas for the spa- 
tial distribution of the order parameter and the magnetic 
field agree with the numerical results to an accuracy of a 
few percent in a wide range of k ^ 1 and B. By introduc- 
ing the effective-mass tensor the theory is extended to in- 
clude anisotropy when the vortices are directed along one 
of the principal axes of the crystal. The results for the lo- 
cal order parameter and the magnetic field are then used 
to calculate the magnetization. The resulting expression 
for the magnetization is in agreement with the London 
model in small fields at k ^ 1, and with the Abrikosov 
approximation at H ~ i?c2- The field dependences of 
the magnetization found by the variational and numerical 
approaches practically coincide. At the same time, the 
difference between the numerical result and the Clem- 
Hao approach is considerable. We also found the field 
dependence of the magnetization in the Wigner-Seitz ap- 
proximation close to Hc2 , where the GL equations can be 
linearized |q|. The calculated magnetization curves are 
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compared with the available experimental data on some 
high-Tc superconductors. 



II. THEORETICAL FORMALISM 

In the Wigner-Seitz approximation both the order pa- 
rameter and the magnetic flux density within the cell 
have axial symmetry. In this case the order parameter 
l^Pl can be presented as /(r) exp(— with radius vec- 
tor r and phase angle if. The free energy density of a 
superconductor can be written as the sum of two contri- 
butions: F — Fern + Fcore- F^rn IS related to the energies 
of magnetic field and supercurrent, and Fcore to the sup- 
pression of IV'I in the vortex core. It is easy to show that 
in the framework of the GL theory in the Wigner-Seitz 
approach Fem and Fcore are given by: 



F — 

^ core — 



27r 



kA, 



cell JQ 



2tt 



kA 



cell JO 



Kr 



rdr, (4) 



rdr, (5) 



where a is the dimensionless vector potential, R and 
Aceii = TT-R^ are the cell radius and area, related to the 
magnetic induction B by Aceii = 2tt / Bk. 
The two GL equations can be written as: 



K r dr V dr 



(7) 



The magnetic field and the vector potential are related 

by 



h = 



1 d{ra) 
r dr 



(8) 



These equations must be supplemented by the boundary 
conditions for the magnetic field and the order parame- 
ter: 



h{R) = he, 



/(O) = 0, f{R) = 0, 



o , . dh 

rf ir)— = -lhi at r ^ 0. 
^ ' dr ' 



(9) 



(10) 



(11) 



Condition (11) follows from Eqs. (7) and (8) |]. The 
system of Eqs. (6)-(ll) is much simpler than the similar 
equations for the hexagonal unit cell. However, even this 



system can be solved only numerically Q . Nevertheless, 
in high and small fields some results can be obtained 
analytically. At small fields, the approximate solution 
to the GL equations in the Wigner-Seitz approach was 
found in Ref. [Q. Here, the spatial variation of the order 
parameter at k 3> 1 can be neglected when calculating 
the magnetic flux density. In this case, h can be found 
analytically from the second GL equation (6) and the 
boundary conditions (9) and (11). This yields the mag- 
netization 10]: 



AtiM{B) = Hci 



1 

2k 



-t- 



1 



Ki{R) 

h{R) ' 2/2 (i?) 



B. (12) 



When H ^ Hd, the radius of the cell is i? ^ 1, and 
Eq. (12) can be expanded in powers of R, yielding: 

- AttM{B) = Hci-^ [\n{2K{H - Hd) + a] , (13) 

with (7 = 1.3456. A similar relationship was obtained in 
the London limit in Ref. |^ for the regular FLL, with 
a = 1.3431 1^ for the triangular lattice. 

At high fields, the magnetization is given by |^: 



M = 



H -H, 



c2 



(14) 



47r/3A(2K2 - 1)' 
where (3a depends only on the symmetry of the FLL 



Pa 



[jPd^i 



(15) 



Here the integrals are taken over the area of the unit cell. 
Let us find the value of (3^ for the circular cell. Near Hc2, 
the magnetic field undergoes only slight spatial variation. 
The vector potential in this case is a « Br/2. The order 
parameter is small at H ^ Hc2, so the first GL equation 
(6) can be linearized. The resulting equation has the 
analytical solution pi| : 



fir) 



s exp 



KrB\ fB 



K KrB 



(16) 



where $ is the Kummer function. Factor s depends on 
the nonlinear term in Eq. (6) but it does not affect the 
value of [3a- Using Eqs. (15) and (16), we find f3A = 
1.1576 for the circular cell. This value is close to [3a — 
1.1596 calculated in Ref. ||3|] for the triangular lattice. 

Thus, in both limits of low and high fields, the magne- 
tization in the Wigner-Seitz approximation is in good 
agreement with that for the regular triangular FLL, 
which has the lowest energy. In the next section, we pro- 
pose a variational model to solve the GL equations in the 
whole field range between Hd and Hc2 at any k > 1/ \/2 
with good accuracy. 
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III. VARIATIONAL PROCEDURE 



Using the trial function (1) multiplied by a variational 
prefactor, /(r) = jaoTj \Jr'^ + Q,, allows us to solve the 
second GL equation (7) analytically within the Wigner- 
Seitz cell: 



h[r) = uW^^Jr^T^) + i;ifo(/oo V^^^+^f), (17) 



^v{B,k) = X 




where u and v can be found from the boundary conditions 
(9) and (11): 



fo 



KlifooP) 



K^v Ki{faoiv)Il{foap) - h{fooiv)Ki{Jocp) 



, (18) 



K^v Ki{fooS,v)h{fooP) - h{foaiv)Ki{focP) ' 



(19) 



and we introduced the notation p = \/B?+Q,. Note 
that, rigorously speaking, the Clem trial function (1) does 
not meet the condition (10) since its derivative can not 
be equal to zero at the cell boundary. However, df/dr is 
small at r = i? and the comparison between the results 
of variational and numerical methods demonstrates good 
accuracy of the approach. 

The values of variational parameters should be found 
by minimization of the total free energy density F — 
Fern + Fcore- Using Eqs. (7) and (4), it is possible to 
obtain the following expression for the magnetic energy 
density: Fem = Bh{0) ||2^. Taking into account Eqs. 
(17)-(19), we get: 



^ p.m. — 



X 

Ko{f^£.v)Il{.focp) + h{fooiv)Ki{foop) 

i^i(/ooe.)/i(/ocp)-/i(/cx.e.)ifi(/oop)' ^ ^ 



The energy related to the spatial variation of the order 
parameter in the vortex core is found from Eq. (5) by a 
straightforward integration [p3|j26| : 



Fcore ~ 



In 



1 



I OC 



Bfl (1 + 



2 2 + BK^i ^{2 + BKil) 



(21) 



The field dependence of the variational parameters is 
calculated numerically by minimization of the total free 
energy F{B, K,£^y, foo) with respect to ^„ and foo- This 
dependence can be approximated by the following ana- 
lytical expressions with an accuracy of about 0.5%: 



f^iB,K)= 1- 




(23) 



where t = 0.985, ^^,0 is the value of at _B = 0. The 
latter can be calculated from the condition dF/d£^y — 
at B = 0: 



(24) 



When K > 1, Eq. (24) has the solution ^^,0 = \/2/k (\/2^ 
in dimensional variables). Eqs. (1), (17)-(19), (23), and 
(24) give the distributions of the order parameter and the 
magnetic field within the Wigner-Seitz cell. 




FIG. 1. The spatial distribution of the dimensionless order 
parameter in the unit cell of the flux-line lattice at difi'erent 
magnetic inductions B = 0.1Hc2, B = 0.5Hc2, B — 0.8Hc2 for 
K — 10. The solid lines correspond to the variational calcu- 
lations in the Wigner-Seitz approximation. The dashed lines 
correspond to the numerical solution of the Ginzburg-Landau 
equations for the triangular lattice (nearest neighbor vortices 
are in the plane of the graph). The distance is measured in 
units of the intervortex spacing d in the triangular lattice. 

The upper critical field is determined as the field at 
which the order parameter in the superconductor be- 
comes equal to zero. As can be seen from Eq. (23), 
one has f = aX H = 0.985k. Thus, the difference be- 
tween the exact Hc2 ~ t and its calculated value is about 
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1.5%. This result is quite natural, since variational pro- 
cedures in general give only approximate solutions to the 
GL equations. Similarly, the Clem value of ff ^i s lightly 
differs from the numerically calculated one ||l^Jl9|] . 

Now we compare the obtained results for the order 
parameter and the magnetic field with the similar dis- 
tributions computed for the triangular lattice by the nu- 
merical method proposed in Ref. ||l6|| . The dependence 
of the order parameter on the distance from the cell cen- 
ter is shown in Fig. 1 for different values of the magnetic 
induction at k = 10. The spatial distribution of the or- 
der parameter in the triangular lattice along the nearest 
neighbor direction is also shown in Fig. 1. The results 
of our variational calculations are close to the numerical 
ones at any values of the magnetic induction. The differ- 
ence does not exceed several percent. Such an accuracy 
of our approach remains in a wide range of k ^ 1 . The 
magnetic flux density as a function of the distance from 
the vortex axis is shown in Fig. 2 a,t k = 10, B — 0.5i?c2- 
There is good agreement between the variational and nu- 
merical results for h. In Fig. 3, the spatial average of 
the order parameter squared, cu = (IV'P), is plotted as 
a function of the magnetic induction at k = 100. The 
comparison with the numerical result shows that the de- 
viation of variationally-calculated to from the exact de- 
pendence does not exceed one percent in small and in- 
termediate fields. Near Hc2 this deviation increases due 
to a small difference between the calculated Hc2 and the 
exact value. 
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FIG. 2. The spatial distribution of the dimensionless mag- 
netic field in the unit cell of flux-line lattice at the magnetic 
induction B = 0.5Hc2 for k = 10. The solid line corresponds 
to the variational calculations in the Wigner-Seitz approxima- 
tion. The dashed line corresponds to the numerical solution. 

Thus, the results of our variational approach agree well 
with the exact numerical solution to the GL equations. 
In the the next section we shall apply this approach to 
the calculation of the magnetization curve. 
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FIG. 3. Averaged order parameter squared lo vs B in di- 
mensionless units for k — 100. The solid line corresponds to 
the variational calculations and the dot line to the numerical 
solution. 

The anisotropy in the GL functional can be taken into 
account by introducing the phenomenological effective- 
mass tensor rrij (j = 1,2,3), where rrij are the effec- 
tive masses in the direction of the principal axes Xj. It 
was shown in Ref. p3[ that the Ginzburg-Landau equa- 
tions can be transformed to isotropic form by a simple 
transformation if k is replaced by k* — npi where 
pi — mi/^TOir7i2m3, (the vortices are directed along the 
xi axis). Thus, in this case, the Wigner-Seitz approxi- 
mation can be used. The Wigner-Seitz cell has an ellip- 
tic shape and can be transformed to the circular cell by 
a scaling transformation: the distance along the X j axis 
should be normalized by mim2'm^ / ^mimj . The magne- 
tization of the anisotropic superconductor can be found 
from the magnetization of the isotropic one by replacing 
the GL parameter k by k*. 



IV. MAGNETIZATION 



The magnetization is defined by the well-known rela- 
tionship: 



A-kM = H - B. 



(25) 



It can be calculated by two equivalent methods if the 
exact solution to the GL equations is known. First, the 
magnetic field H may be calculated by minimization of 
the Gibbs free energy G = F- 2BH: 



IdF 
^~2dB 



(26) 



This derivative was calculated, e.g. in Refs. [| 23|. The 
second approach uses the virial theorem for the flux-line 
lattice, which was proven in Ref. [ p4| , namely the applied 
magnetic fleld H can be found from the local magnetic 
field h and the order parameter / as 
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H = 



2BA 



cell 



J {f-f + 2h')d\ 



(27) 



where the integral is taken over the area of the unit ceU. 
Both methods are equivalent if the exact solutions / and 
h are used Q . 

The variational model gives the spatial distributions 
of the order parameter and the magnetic field within the 
cell, which are close to the exact results. Let us find the 
magnetization by means of both methods. According to 
Eq. (26) the magnetization M is: 



4ttM = 



foa 



, Ko{foo^^)Ii{foop) + /o(/oogt,)i^l(/oop) 
' K^{fUv)h{fooP) - h{foo£.v)K^{.fooP) 

-/l(/ooC.-)^l(/ooP)}-' + 



2 , (2 + 3B<2) 



2 f2 e2 1 f2 
^ J oo'^v r ^ J OO 



/^(2 + 3i?<2) /^(2 + 3B<2) 



2k(2 + B<2)J 2K{2 + BKQy 



}• (28) 



Thus, in the former case the dependence of the magneti- 
zation on the magnetic induction B is given by Eqs. (28) 
and (22)- (24). The dependence of iJ on B is given by 
Eq. (25). Thus, we find the implicit function M{H). 

Within the second approach, the integral (27) can be 
calculated only numerically when / and h are defined by 
Eqs. (1) and (17). Our calculations show that not only 
the values of H found by both methods coincide, but also 
the values of M , which is usually much smaller than H , 
are practically indistinguishable. The difference between 
them is much less than one percent at any induction and 
K > l/\/2. Note that the result for M in the second 
approach is very sensitive to the perturbations of /(r) 
and h{r). For example, if one puts /oo = 1 near Hd and 
minimizes the free energy only with respect to ^„ this 
does not change /(r) and h{r) considerably. However, 
this procedure would lead to an appreciable change of 
M{H) when using Eq. (27), while according to Eq. (26) 
the magnetization practically remains the same. Below, 
we shall use Eq. (28) for the magnetization. 

At low fields the variational parameters (22) and (23) 
may be considered as constants independent of B when 
K ^ 1. In this case, Eq. (28) can be expressed as a 
power series in terms of As a result, it is possible to 
obtain Eq. (12) with H^i given by Eq. (3) at £ = 0.52. 
Thus, in small fields the model reduces to the London 
approximation provided that the variationally-calculated 
value of Hci is used, which is practically indistinguish- 
able from the exact Hd- Actually, the use of the exact 
Hci in small fields is equivalent to taking into account 
the effect of vortex cores. The field dependence of the 



magnetization (28) is shown in Fig. 4 for k = 100. The 
magnetization curves corresponding to the London and 
Abrikosov approximations are also plotted. At low fields, 
the magnetization practically coincides with the results of 
the London approach, the difference between them does 
not exceed 0.5%. AX H ^ Hc2, the behavior of the mag- 
netization is in good agreement with the Abrikosov high- 
field result (14). Although near Hc2 both curves are close 
to each other, the error of our approximation is not so 
small as in low fields because of the slight deviation of 
the calculated Hc2 from the exact value. For example, 
at H = 0.8Hc2, K — 100 the error of our variational pro- 
cedure is about 5%. In order to improve the accuracy 
near Hc2, one may put the constant t in Eq. (23) to be 
equal to 1. As a result, the difference between the varia- 
tional and the Abrikosov M{H) curves decreases in the 
vicinity of i?c2, whereas at low and intermediate fields 
the magnetization does not change. 

In the inset of Fig. 4, we compare the calculated de- 
pendence —A'kM{H) with that found in Ref. where 
the Clem-Hao approach of superposition of vortex fields 
p3[ was used. It is clearly seen that this approach is valid 
only at small fields, and its use leads to an appreciable 
error in the magnetization even in the intermediate field 
range. An additional error in M arises in small fields 
due to the approximate replacement of the lattice sums 
by integrals psf in calculating the magnetic free energy; 
for more details see Ref. [BOl. 
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FIG. 4. Calculated —AtiM{H) using the variational 
method for k = 100 (solid line). Also shown are the Lon- 
don dependence —A-nM{H) (13) with the exact value of Hci 
(3) (dashed line) and the Abrikosov high- field result (16) (dot 
line). The inset refers to the magnetization found within the 
framework of the variational model (solid line) and the result 
of calculations teo| for the magnetization in the Clem-Hao 
approximation [p3| (dot line). The difference between these 
curves arises due to the superposition of vortex fields used in 
the Clem-Hao model |^^. Dimensionless variables are used. 

Now we compare the magnetization found in the frame- 
work of the variational and numerical methods. At low 
and high fields the exact dependence M{H) coincides 
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with the results of the London (at k ^ 1) and the 
Abrikosov approximations, respectively. As we found 
above, the results of the variational approach are in 
agreement with these approximations. In the interme- 
diate field range, where the London and the Abrikosov 
approaches are not applicable, the difference between the 
values of the magnetization calculated by numerical and 
variational methods is not bigger than 1% in a wide range 
of K ^ 1 values. Thus, our results for the magnetization 
appear to be a good approximation to the exact numer- 
ical solution of the GL equations at k 1. 

Next we discuss the case of small k values. In Fig. 5, 
the field dependences of the magnetization are plotted for 
several small k values. The solid and dotted lines corre- 
spond to the variational and numerical calculations, re- 
spectively. The agreement between these results is good 
at K ^ 1. At smaller k values the variational and the 
exact numerical results differ near the lower critical field. 
In this case, the intervortex distance is of the order of the 
coherence length almost in the entire field range, and the 
variational approach based on an appropriate trial func- 
tion for the order parameter in the circular Wigner-Seitz 
cell may lead to some deviation from the exact solution. 

In the inset of Fig. 5 we compare the magnetization cal- 
culated by means of our variational procedure and by the 
Clem-Hao model at k = 1. It is seen that also for small 
(even as for large k, see above) the calculation method 
proposed in Ref. |23| leads to an inaccurate magnetiza- 
tion. 



analyzed and reduced to the dimensionless form based on 
the Clem-Hao formulas with non-selfconsistent field de- 
pendences of the variational parameters . The result- 
ing magnetization curve is close to the Abrikosov high- 
field result and to our variational dependence in the in- 
termediate field range. The n values obtained in Refs. 
||,|9| were k = 70 for YBa2Cu408 and k = 80 for 
Ndi.85Ceo.i5Cu04-5. 



0,010 




0,000 



FIG. 6. The field dependence of magnetization in dimen- 
sionless units. The solid and dotted lines show the theoretical 
variational dependences at k = 70 and «: = 80. The circles 
and triangles give the experimental data for YBa2Cu408 |2^ ] 
and Ndi.85Ceo.i5Cu04_5 respectively. The inset shows 
the magnetization calculated by our variational model and 
plotted versus the logarithm of the applied field H at k — 100. 




FIG. 5. Magnetization curves calculated using the vari- 
ational approach (solid lines) and exact numerical method 
(dotted lines) at various small k values. The inset compares 
the magnetization obtained by our variational model (solid 
line) with the Clem-Hao approximation (dashed line) at 
K, — 1. Dimensionless variables are used. 

Our formulas for the magnetization may be used for 
the analysis of experimental data. In Fig. 6 the calcu- 
lated magnetization curves are compared with the mea- 
sured magnetization of YBa2Cu408 polycrystals and 
Ndi.g5Ceo.i5Cu04_5 single crystals |29j. In these papers, 
the magnetization curves at different temperatures were 



In Fig. 6 we compare these experimental curves with 
our variational result. The circles and triangles in Fig. 6 
correspond to YBa2Cu408 and Ndi.85Ceo.i5Cu04_5, re- 
spectively. The solid and dotted lines correspond to the 
theoretical curves at k = 70 and k = 80, respectively. It 
is clearly seen that good agreement exists between the- 
ory and experiment. In the inset of Fig. 6 the theoretical 
(variational) magnetization is plotted as a function of 
InH at K — 100. This dependence is nearly linear in 
a wide range of intermediate fields. A similar behavior 
of the magnetization was observed in numerous experi- 
ments; see, for example, Ref. |20| . 

V. CONCLUSIONS 

We proposed an approximate method to solve the 
Ginzburg-Landau equations for the regular flux-line lat- 
tice at any values of the magnetic induction and the 
Ginzburg-Landau parameter k > 1/V2. The Wigner- 
Seitz approximation is used, and the hexagonal unit cell 
of the vortex lattice is replaced by a circle with the same 
area. Our model is based on Clem's trial function for 
the order parameter. The use of this function allows us 
to find the magnetic flux density self-consistently from 
the second Ginzburg-Landau equation. The comparison 
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between the variational results and the results of exact 
numerical solution of the Ginzburg-Landau equations re- 
veals good accuracy of our approach: the difference be- 
tween the spatial distributions of the order parameter 
and the magnetic field does not exceed several percent. 
Such accuracy remains in a wide range of values of the 
magnetic induction and 1. The method is applied to 
the calculation of the field dependence of the reversible 
magnetization. An analytical expression for the mag- 
netization is proposed. At low fields, the obtained de- 
pendence agrees with the predictions of London theory 
at K ^ 1. At high fields, it is in good agreement with 
the Abrikosov result. It is shown that the values of the 
magnetization calculated within the framework of our 
variational model and of the numerical method of solu- 
tion of the Ginzburg-Landau equations are practically 
indistinguishable, especially in small and intermediate 
magnetic fields at k > 1. Our model yields the lim- 
its of the Clem-Hao model for the magnetization. The 
presented analytical formulas for the magnetization may 
be used to analyze experimental data. As an illustra- 
tion we compared the experimental and calculated mag- 
netization curves for different high-Tc superconductors 
(YBa2Cu408 and Ndi.85Ceo.i5Cu04_i) and found good 
agreement between theory and experiment. 
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